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Abstract 

A  laser  deposition  process  involves  the  supply  of  metallic  powders  into  a  laser-heated  spot  where 
the  powder  is  melted  and  forms  a  melt  puddle  which  quickly  solidifies  into  a  bead.  In  order  to 
design  an  effective  system,  the  laser  beam,  the  powder  beam,  and  their  interactions  need  to  be 
fully  understood.  In  this  paper,  the  laser-material  interaction  within  the  melt  pool  is  reported 
using  a  multi-scale  model;  a  macroscopic  model  to  model  mass,  heat,  and  momentum  transfer. 
Experiments  were  also  conducted  to  validate  the  simulation  model. 

Introduction 

The  direct  laser  deposition  (DLD)  process  is  an  extension  of  the  laser  cladding  process  for  rapid 
prototyping  of  fully  dense  metal  components.  It  offers  the  ability  to  make  a  metal  component 
directly  from  CAD  drawings.  It  involves  the  supply  of  metallic  powders  into  a  laser-heated  spot 
where  the  powder  is  melted  and  forms  a  melt  puddle  which  quickly  solidifies  into  a  bead  as 
shown  in  Figure  1.  The  thermal  behavior  of  the  melt  pool  in  which  powder  is  injected  is  critical 
for  part  quality.  This  paper  uses  simulation  and  experiments  to  investigate  the  effect  of  the 
process  parameters:  laser  power,  powder  mass  flow,  and  scanning  speed  on  the  melt  pool  thermal 
behavior.  During  the  laser  deposition  process,  several  defects,  such  as  porosity  and  cracks, 
should  be  paid  attention  to.  Cracks  initiate  corrosion  fracture  and  reduce  fatigue  strength  of  the 
deposited  parts.  Cracks  are  caused  by  the  residual  stresses  created  by  the  high  thermal  gradient 
built  up  during  the  cooling  stage.  Residual 
stresses  can  be  reduced  by  a  reduction  of 
the  cooling  rate.  This  can  be  achieved  by 
preheating  the  substrate.  Moreover,  the 
preheating  causes  a  better  absorption  of  the 
laser  beam  and  further  it  is  possible  to  melt 
more  powder  in  the  larger  melting  pool  and 
enhance  the  bonding.  Usually  the 
preheating  is  performed  in  such  a  manner 
that  a  very  small  melting  of  the  substrate 
material  occurs  without  powder  injection. 

Thermal  analysis  of  the  laser  deposition 
process  is  very  important  for  optimization 
of  the  process.  If  the  substrate  surface 
temperatures  remain  too  low,  wetted 
deposit  material  is  limited.  In  that  case. 


Fig  1  Laser-powder  interaction  at  melt  pool  in  a 
direct  laser  deposition  (DLD)  process. 
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irregularly  shaped  traeks  with  a  lot  of  cracks,  porosity  and  a  poor  bonding,  are  produced. 
However,  if  an  extremely  high  temperature  of  the  substrate  surface  is  reached,  severe  melting  of 
the  substrate  occurs.  The  high  degree  of  dilution  can  deteriorate  the  clad  properties.  A 
comprehensive  numerical  model  has  been  developed  that  allows  the  prediction  of  temperature 
distribution  and  melt  pool  dynamics.  This  model  simulates  the  coaxial  laser  deposition  process 
with  powder  injection,  and  considers  most  of  the  associated  phenomena,  such  as  melting, 
solidification,  evaporation,  evolution  of  the  free  surface,  and  powder  injection.  Input  parameters 
for  this  model  are  laser  machining  parameters  and  properties  of  the  laser  beam,  as  well  as 
material  properties  and  laser  beam  absorption.  To  get  more  accurate  predictions,  a  finer  grid 
needs  to  be  used.  This,  together  with  the  iterative  nature  of  the  numerical  algorithms,  causes  the 
model  computationally  not  to  be  very  efficient.  In  this  paper  an  analytical  model  is  applied  to  the 
preheating  process  (without  powder  injection)  to  increase  the  computational  efficiency,  while  the 
actual  deposition  process  with  powder  injection  still  uses  the  numerical  model.  The  outputs,  and 
the  temperature  distribution  of  the  substrate,  are  used  as  the  initial  conditions  of  the  numerical 
model. 

In  this  study,  a  coaxial  diode  laser  deposition  system,  LAMP  (developed  at  UMR),  is  considered 
for  simulations  and  experiments.  The  blown  powder  method  is  used  to  deliver  powder.  As  the 
primary  laser  in  the  LAMP  system,  the  diode  laser  is  used.  Material  of  both  powder  and 
substrates  is  Ti-6A1-4V,  which  is  widely  used  in  the  aerospace  industry. 

Melt  Pool  Governing  Equations 

Figure  2  shows  a  schematic  diagram  of  the  calculation  domain,  including  the  substrate,  melt 
pool,  remelted  zone,  deposited  layer  and  part  of  the  gas  region.  In  the  laser  deposition  process, 
melting  and  solidification  cause  the  phase  transformation  at  the  solid/liquid  interface.  A  mushy 
zone  containing  solid  and  liquid  is  formed.  In  this  study  the  continuum  model  [BennonSVa, 
BennonSVa]  is  adopted  to  derive  the  governing  equations. 


Fig.  2.  Schematic  of  the  calculation  domain. 


The  assumptions  for  the  system  of  governing  equations  include:  (1)  the  fluid  flow  in  the  melt 
pool  is  a  Newtonian,  incompressible,  laminar  flow;  (2)  the  solid  and  liquid  phases  in  the  mushy 
zone  are  in  local  thermal  equilibrium;  (3)  the  solid  phase  is  rigid;  and  (4)  isotropic  permeability 
exists.  For  the  system  of  interest,  the  conservation  equations  for  mass,  momentum  and  energy  are 
summarized  as  follows: 
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Continuity 


^  +  V.(pf)  =  0 

at 


Momentum 


^{pV)  +  V  ■{pVV)  =  V  ■{^,-^VV)-Vp-^-^{V -t)  +  pg  +  S 
dt  p^  K  p^ 


Energy 

+  V.{pVh)  =  V-{kVT)-V- {p{h,  -  hW  -  V, ) 


(1) 


(2) 


(3) 


In  equations  (l)-(3),  the  eontinuum  density,  speeific  heat,  thermal  conductivity,  vector  velocity, 
and  enthalpy  are  defined  as  follows: 

p  =  s,p,+giPi  c  =  /a+/c,  *  =  gA+gA 

V  =  ff.  +  f,V,  k  =  f,k,+fA  (4) 


The  liquid  fraction  temperature  relationship  is  given  by: 


gi  = 


0  if  r  <  r 
T-T 

'  if  r  <  r  <  E 


T,-T 

I  s 


1  if  r  >  t; 


The  other  volume  and  mass  fractions  can  be  obtained  by: 


(5) 


f,= 


Ml 

P 


_  ?>sPs 

p 


g,+g/=l  fs=^-fl 


The  phase  enthalpy  for  the  solid  and  the  liquid  can  be  expressed  as: 

h^=[c^(T)dT  h,  =  \\^{T)dT  +  \lc,(T)dT  + 


(6) 


(V) 


where  Lm  is  the  latent  heat  of  melting. 

Permeability,  K,  is  assumed  to  vary  with  liquid  volume  fraction  according  to  the  Carman  - 
Kozeny  equation  [Carman37]  derived  from  Darcy’s  law: 
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(8) 


K  = 


g 


3 

/ 


where  the  parameter  C  is  a  eonstant  depending  on  the  morphology  and  size  of  the  dendrites  in 
the  mushy  zone.  The  S  is  a  source  term  that  is  defined  below. 

Solid/liquid  interface 

The  solid/liquid  interface  is  implicitly  tracked  by  the  continuum  model  [Bennon87a, 
Bennon87a].  In  the  solid  phase  region  and  liquid  phase  region,  the  third  term  on  the  right-hand 
side  of  Eq.  (2)  vanishes.  This  is  because  in  the  solid  phase  region  f  =  =0  and  in  liquid  phase 

region  K^oo  since  gi  =  1 .  So  this  term  is  only  valid  in  the  mushy  zone. 

Free  surface 

The  liquid/vapor  interface,  or  the  free  surface  of  the  melt  pool,  is  very  complex  due  to  surface 
tension,  thermocapillary  force,  and  impaction  of  the  powder  injection.  In  this  study,  the  Volume- 
Of-Fluid  (VOF)  method  [Hirt8I]  is  employed  to  track  the  evolution  of  the  moving  free  surface  of 
the  melt  pool.  The  melt  pool  configuration  is  defined  in  terms  of  a  volume  of  fluid  function, 
F(x,y,t),  which  represents  the  volume  of  fluid  per  unit  volume  and  satisfies  the  conservation 
equation: 

O  j-r 

—  +  {V-y)F  =  0  (9) 

dt 


Source  term 

The  source  term,  S,  in  the  momentum  equation  is  contributed  by  the  interface  forces  acting  on 
the  free  surface,  such  as  surface  tension,  etc.  In  this  study,  the  continuum  surface  force  (CSF) 
model  [Brackbill92]  is  used  to  reformulate  the  surface  force.  In  its  standard  form,  surface  tension 
is  formulated  as  [Brackbill92]: 

=  +  sY  (10) 


where  F^{xY)  is  the  net  surface  force  at  a  point  on  a  interface  S.  n  is  a  unit  normal  to  S  at  the 
point  ,  which  is  given  by: 

.  n 

n  =  — 

n 

where  n  is  the  surface  normal  vector  and  can  be  computed  from  the  gradient  of  the  VOF 
function: 

n=VF  (12) 

is  the  gradient  along  a  direction  tangential  to  the  interface,  which  is  defined  as: 
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=  V  -  Vjy  =  V  -  h{h  ■  V) 


(13) 


Y  and  k  represent  surface  tension  coefficient  and  curvature,  respectively,  k  is  given  in  [ 
Prakash89]; 


.V.  1 

K  -  -(V  -n)-  — 

n 

rl 

(14) 


By  using  the  CSF  model,  the  surface  force  is  reformulated  into  a  volume  force  as  follows: 


F^n 


F 

F 


F,VF 


F 

F 


(15) 


where  (f)  is  the  averaged  F  value  across  the  free  surface.  Thus,  the  source  term  S  in  equation 
(2)  is  formulated  as: 


5  = 


(ny/f  +  V^y)|VF| 


(16) 


Boundary  Conditions 

The  boundary  conditions  at  the  free  surface  satisfy  the  following  equation: 


dn  71 


-hXT-TJ-sa{r-T:)-mL^ 


(17) 


where  terms  on  the  right-hand  side  are  laser  irradiation,  convective  heat  loss,  radiation  heat  loss 
and  evaporation  heat  loss,  respectively.  Piaser  is  the  power  of  laser  beam,  Patten  the  power 
attenuated  by  the  powder  cloud,  R  is  the  laser  beam  radius,  p  the  laser  absorption  coefficient. 
Patten  is  Calculated  according  to  Frenk’s  et  al.  model  [Frenk97]  with  a  minor  modification: 


P  =P 

atten  laser 


( 


1-exp 


(18) 


where  m  denotes  the  powder  mass  flow  rate,  1  is  the  stand-off  distance  from  the  nozzle  exit  to  the 
substrate,  p  is  powder  density,  rp  is  the  radius  of  the  powder  particle,  Djet  is  the  diameter  of  the 
powder  jet,  Vp  is  the  powder  injection  velocity,  and  Qext  is  the  extinction  coefficient.  It  is 
assumed  that  the  extinction  cross  section  is  close  to  the  actual  geometrical  cross  section,  and  Qext 
takes  a  value  of  unity.  In  the  evaporation  term,  is  the  evaporation  mass  flux  and  Ly  is  the 
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latent  heat  of  evaporation.  Aoeording  to  Choi  et  al.’s  overall  evaporation  model  [Choi87],  is 
of  the  form; 

logm^  =A  +  6.1210-^^P^-0.51ogr  (19) 

where  A  is  a  eonstant  dependent  on  the  material. 

The  boundary  eonditions  at  the  bottom,  left  and  right  wall  satisfy  the  following  equations: 


dT 

k-^-h  (T-T  ) 

^  c  V  00  / 

on 

(20) 

u  =  0  V  =  0 

(21) 

Powder  injection 

Powder  particles  that  inject  onto  the  top  surface  can  be  classified  into  three  categories.  1)  Those 
powder  particles  that  have  not  been  melted  during  their  passage  and  hit  the  solid  part  of  the 
substrate  will  deflect  and  lose;  2)  Those  powder  particles  that  have  been  melted  before  they 
arrive  on  the  substrate  and  impact  the  solid  part  of  the  substrate  will  stick  to  the  surface  of  the 
substrate;  and  3)  Those  powder  particles  that  fall  into  the  melt  pool  including  melted  and 
unmelted  will  merge  and  mix  with  the  molten  liquid  in  the  melt  pool  [Han04].  In  this  study,  the 
powder  particles  belonging  to  the  last  two  cases  were  modeled. 


Numerical  simulation 

Simulation  is  performed  based  on  the  capability  of  the  experimental  facilities  to  compare  the 
simulation  results  with  the  experimental  measurements.  A  continuous  wave  diode  laser  with  an 
808  nm  wavelength  is  considered  as  the  energy  source.  The  laser  intensity  distribution  is 
uniform.  For  substrates,  Ti-6A1-V4  plates  with  a  thickness  of  0.25  inch  are  selected.  Ti-6A1-V4 
powder  particles  with  a  diameter  from  40  to  140  qm  are  used  as  deposit  material.  The  laser 
absorption  coefficient  is  measured  by  Sparks  et  al.  [Sparks06].  The  material  properties  and  the 
main  process  parameters  are  shown  in  Table  1.  Figure  3  shows  typical  simulation  results  for 
normal  deposition  and  lack  of  fusion. 


Fig.  3.  Typical  results  for  simulation  of  normal  deposition  (left)  and  lack  of  fusion  (right). 
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Experimental  Study 

The  experiments  were  performed  on  the  LAMP  system  which  consists  of  a  diode  laser,  powder 
delivery  unit,  5-axis  CNC  machine,  and  monitoring  subsystem.  The  laser  system  used  in  the 
study  was  Nuvonyx  (Nuvonyx  Inc.)  ISL-IOOOM  Laser  Diode  System  which  combines  state-of- 
the-art  micro-optics  with  laser  diodes  to  produce  the  only  single  wavelength  fiber  coupled  direct 
diode  laser  at  power  levels  up  to  1000  watts  CW.  The  laser  emits  at  808  nm  and  operates  in  the 
continuous  wave  (CW)  mode.  The  substrates  have  dimensions  of  2. 5x2. 5x0.4  in.  The  Ti-6A1-4V 
samples  were  irradiated  using  a  laser  beam  with  a  beam  spot  diameter  of  2.5  mm. 


Table  1.  Material  properties  for  Ti-6A1-4V  and  main  process  parameters 


Nonmenclature 

Symbol 

Value  (unit) 

Melting  temperature 

T™ 

1900.0K 

Liquidus  temperature 

Ti 

1923.0K 

Solidus  temperature 

Ts 

1877.0K 

Evaporation  temperature 

Tv 

3533.0K 

Solid  specific  heat  at  constant 
pressure  [Kelly04] 

Cps 

1 

483.04  + 0.2 15r  T  <126^ 

J !  kg  K 

412.7  +  0.18017’  1268<7’<1923 

Liquid  specific  heat  at  constant 

Cpl 

831.0  J/kgK 

pressure  [Mills02] 

Thermal  conductivity  [ 

k 

1.2595  +  0.01577’  7’<12687: 

Kelly04] 

< 

3.5127+0.01277’  1268<7’<1923  WImK 

-12.752  +  0.0247’  7’ >1923 

Solid  density  [MillsO] 

Ps 

4420-0.154  (T-298K) 

Liquid  density  [MillsO] 

Pi 

3920-0.68  (T-  1923  K) 

Latent  heat  of  fusion  [MillsO] 

Lm 

2.86  X  10^  J/kg 

Latent  heat  of  evaporation 

Lv 

9.83  X  lO'^J/kg 

Dynamic  viscosity 

P 

3.25  X  10'^  N/m  s  (1923K)  3.03  x  10'^  (1973K) 
2.66  X  10'^  (2073K)  2.36  x  10'^  (2173K) 

Radiation  emissivity  [Lips  05] 

s 

0.1536  +  1.8377  x  lO’'*  (T  -  300.0  K) 

Laser  absorption  coefficient  [ 

h 

0.4 

Sparks04] 

Powder  particle  diameter 

Dp 

40-140  pm 

Shielding  gas  pressure 

Pg 

5  psi 

Ambient  temperature 

Too 

300K 

Convective  coefficient 

he 

10  wWk 

In  order  to  validate  the  model  predictions,  single  path  deposition  experiments  are  conducted.  The 
comparisons  between  model  predictions  and  experimental  results  are  conducted  in  terms  of  melt 
pool  peak  temperature.  The  melt  pool  peak  temperature  is  calibrated  through  the  dual¬ 
wavelength  non-contact  temperature  sensor,  which  can  effectively  decrease  the  disturbance  from 
the  powder  and  other  dusts. 

Comparisons  and  discussions 

Fig.  4  and  Fig.  5  show  the  comparisons  between  experimental  measurements  and  model 
predictions.  Fig.  4  shows  the  effects  of  laser  power  on  melt  pool  peak  temperature.  It  can  be  seen 
that  an  increase  in  the  laser  power  will  increase  the  melt  pool  temperature.  This  is  easy  to 
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understand.  As  the  laser  power  inereases,  more  power  is  available  for  melting  the  substrate.  Fig. 
4  shows  the  effects  of  laser  scanning  speed  on  the  melt  pool  peak  temperature.  An  increase  in  the 
laser  scanning  speed  will  decrease  the  melt  pool  peak  temperature.  This  is  because,  as  scanning 
speed  decreases,  the  laser  material  interaction  time  is  extended. 

From  Figs.  4-5,  one  can  see  that  the  general  trend  between  experimental  measurements  and 
model  predictions  is  consistent.  At  a  different  power  intensity  level,  there  is  a  different  error 
from  10  K  (about  0.5%)  to  121K  (about  5%).  It  can  be  seen  that  at  a  higher  power  intensity  level, 
there  is  a  bigger  error  between  measurements  and  predictions.  This  is  because  the  numerical 
model  is  two-dimensional.  It  does  not  consider  the  heat  and  mass  transfer  in  the  third  direction. 
At  a  higher  power  intensity  level,  heat  and  mass  transfer  in  the  third  dimension  are  more 
significant. 

The  errors  between  experimental  measurements  and  model  predictions  are  analyzed  to  mainly 
come  from  the  following  aspects:  (1)  The  two-dimensional  nature  of  the  numerical  model;  (2) 
The  thermo-physical  properties  taken  for  the  analytical  model;  (3)  The  uncertainties  of  the 
material  properties  and  the  appropriateness  of  the  sub-models  for  the  numerical  model;  (4) 
Boundary  conditions.  Adiabatic  boundary  conditions  are  assumed  in  the  analytical  model  and  the 
numerical  model  for  the  bottom  surface  and  side  surfaces.  Measurements  also  have  been  taken  to 
achieve  such  boundary  conditions  in  experiments.  But  it  is  hard  to  get  absolute  adiabatic 
boundary  conditions. 


Laser  Power,  W 


Fig.  4  Comparisons  between  experimental  measurements  and  model 
predictions  at  a  constant  powder  mass  flow  rate  of  5g/min  and  a  constant 
laser  scanning  speed  of  20  inch/min.  (For  actual  deposition) 
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Laser  Scanning  Speed,  inch/m  in. 

Fig.  5  Comparisons  between  experimental  measurements  and  model 
predictions  at  a  constant  powder  mass  flow  rate  of  5g/min  and  a 
constant  laser  power  of  700  W.  (For  actual  deposition) 

The  simulation  model  was  also  used  to  prediet  the  laek  of  fusion  to  deposit  Ti64.  Twenty-eight 
Ti64  samples  were  built  to  validate  the  simulation  model.  The  idea  was  to  use  the  simulation 
model  to  predict  the  lack  of  fusion,  and  the  corresponding  laser  power  is  used  as  nominal  power. 
The  depositions  are  then  based  on  Nominal,  Nominal  +-  10%,  Nominal  +-  20%,  and  Nominal  +- 
30%.  The  nominal  is  about  600  W  for  the  direct  diode  laser.  The  nominal  power  is  then  set  at 
600w  with  a  travel  speed  (scanning  speed)  at  lOipm  and  powder  flow  rate  ~  5g/min.  The 
experiment  was  performed  at  -30%  (420w),  -20%  (480W),  -10%  (540w),  nominal  (600w), 
10%(660w),  20%(720w),  30%  (780w).  If  the  predicted  model  for  lack  of  fusion  is  accurate, 
about  50%  of  the  parts  will  be  observed  as  lack  of  fusion.  Seven  samples  of  thick  walls  are 
shown  in  Figure  6.  It  can  be  observed  that  most  samples  with  lower  power  than  the  nominal 
showed  obvious  porosities  on  the  sample  surface. 


(a)  Sample  produced  at  -30%  laser  power  (b)  Sample  produced  at  -20%  laser  power 
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(c)  Sample  produced  at  -10%  laser  power  (d)  Sample  produced  at  nomiual  laser  power 


(e)  Sample  produced  at  +10%  laser  power  (f)  Sample  produced  at  +20%  laser  power 


(g)  Sample  produced  at  +30%  laser  power 


Fig  6.  A  thick  wall  Ti64  part  produced  at  various  laser  powers  (Nomiual,  Nomiual  +- 10%,  Nomiual  +-  20%, 
aud  Nomiual  +-  30%)  iu  the  UMR  LAMP  lah.  Lack  of  fusiou  cau  he  observed  from  some  surfaces  of  the 
samples  with  lower  power. 


Conclusions 

The  thermal  behavior  of  melt  pool  is  the  most  critical  phenomenon  in  laser  deposition.  An 
analytical  model  for  thermal  analysis  of  temperature  rise  due  to  a  moving  heat  source  is 
combined  into  a  comprehensive  heat  transfer  and  fluid  flow  numerical  model  for  the  laser 
deposition  process.  The  analytical  model  is  used  for  the  preheating  process  before  the  actual 
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laser  deposition  with  powder  injeetion.  And  the  numerieal  model  is  used  for  the  aetual  laser 
deposition  proeess.  Thus  the  outputs  of  the  analytieal  model  are  used  as  the  inputs  of  the 
numerieal  model.  Experiments  have  been  eondueted  to  validate  the  model  predietions.  A 
consistent  general  trend  is  found  between  experimental  measurements  and  the  model  predictions. 
The  sources  of  the  errors  have  also  been  analyzed. 
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